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Abstract 

The quantum analogy of the usual spinodal decomposition is explored for the re- 
cently achieved binary alkali Bose-Einstein condensates mixture. We conclude that 
an analogy is possible within the formulation of coupled non-linear Schrodinger equa- 
tions, and find that the quantum spinodal decomposition consists of two stages. The 
non-equilibrium stage I is dominated by the fastest growth mode, associated with a 
characteristic length. Expressions for both time and length scales in the stage I are 
obtained. The state II is a relaxation process of approaching equilibrium, dominated 
by the slowest mode. We propose that for this state the slow evolution towards the 
phase segregation is due to the Josephson effect between different domains of same 
condensate, and its time scale is estimated. 
PACS#: 03.75.Fi; 64.75.+g 
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I. INTRODUCTION 



Spinodal decomposition in a binary-solution system is a typical example of phase ordering dy- 
namics: the growth of order through domain coarsening when a system is quenched from the 
homogeneous phase into a broken-symmetry phase KJ. Systems quenched from a disordered phase 
into an ordered phase do not order instantaneously. Instead, different length scales set in as the 
domains form and grows with time, and different broken symmetry phases compete to select the 
equilibrium state. In the dynamical equation, the Cahn-Hilliard equation, for the spinodal decom- 
position, the diffusion constant of the material plays a decisive role and determines both the length 
and time scales. It appears that a system without the diffusion constant may not show all the 
main features of the spinodal decomposition. With the recent realization of the binary Alkali Bose- 
Einstein condensates (BEC's) [[| we demonstrate in the present article that indeed it is possible to 
have a spinodal decomposition without the diffusion constant and that this may have been realized 
experimentally. The binary alkali BEC mixtures provide new systems to non-equilibrium phenom- 
ena in new parameter regimes. Its mathematical description is simple enough that a theoretical 
description of the whole process is feasible. A direct comparison between theoretical calculations 
and experimental observations can be made. Here, we shall study the gross features in the dynam- 
ical evolution process, starting from the homogeneous unstable state. To differentiate the present 
situation from the usual ones, we shall call the present one the quantum spinodal decomposition, 
and the previous ones the classical spinodal decomposition. 

In a classical spinodal decomposition process, the particle number for each specie is conserved 
separately. The process can be classified into two stages, the initial stage of fast growth from the 
homogeneous unstable state and the late stage towards to the equilibrium of true ordering. The 
initial stage is dominated by the fastest growth mode determined by the dynamical equation. It 
is a highly non-equilibrium process. There is a length scale associated with this time scale, which 
gives the characteristic domain size after this initial stage. In the late stage the domains grow and 
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merge, and their numbers becomes smaller and smaller. This is a relaxation process towards the 
equilibrium, dominated by the slowest time scale. For an infinitely large system, the time scale may 
be infinity. Hence the system may never achieve true equilibrium. We shall show below that the 
quantum spinodal decomposition occurring in the binary BEC mixtures shares all the features of 
the classical spinodal decomposition. The main difference is that the size of binary BEC mixtures 
is finite. Therefore it is possible to achieve true equilibrium within a finite time scale. To further 
differentiate the present process from the classical one, we shall call the initial stage in the dynamical 
evolution of the binary BEC mixtures the stage I, and call the second stage towards equilibrium 
the stage II. 

In the following, we first formulate the problem in Sec. II. The coupled nonlinear Schrodinger 
equations and the parameter regime to realize the quantum spinodal decomposition will be specified. 
In Sec. Ill the stage I is studied in detail. The fastest growth mode is explicitly given. Interspersed 
regions or domains of coexisting condensate 1 and 2 are formed at this stage, characterized by a 
length scale. Stage II is studied in Sec. IV. We first demonstrate that the scenarios in classical 
spinodal decompositions are not possible here. We propose and analyze a new mechanism for the 
approach towards equilibrium: the Josephson effect between different domains of same condensate. 
This is in accordance with the present quantum spinodal decomposition concept. A comparison 
to experimental situations is discussed in Sec. V. We conclude there that the quantum spinodal 
decomposition can be realized. We summarize in Sec. VI. 



II. FORMULATION OF PROBLEM 



We start from the Hamiltonian formulation of a binary BEC mixture at zero temperature: 
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Here ipj, rrij, Uj with j = 1, 2 are the effective wave function, the mass, and the trapping potential 
of the jth condensate. The interaction between the jth condensate atoms is specified by Gjj, 
and that between 1 and 2 by Gi 2 . In the present paper all G's will be taken to be positive. 
The corresponding time dependent equations of motion are the well-known non-linear Schrodinger 
equations ||, obtained here by minimization of the action, S = J dt{J2j=i,2 ^j^m^Pj ~ -^0' 
d h 2 

ih—iPi(x,t) = -- — V 2 ViOM) + (Ui(x) - fii)iPi(x,t) + Gu\i)i(x,t)\ 2 i)i(x,t) 
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The Lagrangian multipliers, the chemical potentials /ii and fi 2 , are fixed by the relations 
/ d 3 x\ipj(x, t)\ 2 = Nj ,j = 1,2, with Nj the number of the jth condensate atoms. Eq. (2) and 
(3) are mean field equations, since we treat the effective wave functions as c numbers, correspond- 
ing to the Hartree-Fock-Bogoliubov approximation. 

Experimentally, the trapping potentials {Uj} are simple harmonic in nature. For the sake of 
simplicity and to illustrate the physics we shall consider a square well trapping potential Uj = U: 
zero inside and large (infinite) outside, unless otherwise explicitly specified. We will come back to 
this question in the discussion of experimental feasibility. 

We shall consider the strong mutual repulsive regime 



Gi 2 > y GuG 22 . (4) 
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In this regime the equilibrium state for two Bose-Einstein condensates is a spatial segregation of 
two condensates, where two phases, the weakly and strong segregated phases, characterized by the 
healing length and the penetration depth, have been predicted We shall use Eq. (2) and (3) 
under the condition (4) to study a highly non-linear dynamical process: The two condensates are 
initially in a homogeneously mixed state, then eventually approach to the phase segregated state. 
In the same mean field manner as in Ref. we find that this dynamic process can be classified 
into two main stages: The initially high non-equilibrium dynamical growth in stage I, where the 
dynamics is governed by the fastest growth mode, and the stage II of approaching to equilibrium 
where the dynamics is governed by the slowest mode. The stage II is the typical of a relaxation 
process near equilibrium. However, we shall show again it is governed by a quantum effect, namely, 
the Josephson effect. 



III. STAGE I: FASTEST GROWTH MODE 

With the square well trapping potential specified in Sec. II, the coupled non-linear Schrodinger 
equations have an obvious homogeneous solution: Inside the trap the condensate densities \ipj\ 2 = 
Pj , pjQ = Nj/V , with V the volume of the square well potential trap, and the chemical potentials 
Pi = Gupio + G12P20 and \ii = G22P20 + Ci2Pio- This is the initial condition of the present problem. 
It is known that for a large enough mutual repulsive interaction, that is, if G12 > , this 

initial state is not the ground state. Rich physics has been displayed by various theoretical studies 
|5|||. Among many unstable and growing modes in this parameter regime, we will find the fastest 
growth mode in the initial stage of the process towards to the equilibrium state. 

To look for the fastest growth mode out of the homogeneous state, we start with small fluctu- 
ations from the homogeneous state. This is consistent with the usual stability analysis ||. Our 
approach here is to emphasize the connection with the physics of the classical spinodal decomposi- 
tion and the role played by the Josephson relationships. Define 
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and define the density fluctuations Spj = pj — pj and the phase fluctuations Oj, and assume they 
are small: \Spj\/pj , << 1- The definition of the phase fluctuations here has made use of the 
implicit assumption that there is no average current. Then, to the linear order, we have from Eq. 
(2,3), for condensate 1, 

n 
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In terms of hydrodynamics, Eq.(7) is the continuity equation, and Eq.(8) is the Bernoulli equa- 
tion, with the first term at the right hand side as the so-called 'quantum pressure'. Similarly, for 
condensate 2 we have 
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Eliminating the phase variables from Eqs. (7-10), we have 
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We look for the solution of the form 
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with A, B constants. Eq. (11) then becomes 
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In order to have a non-zero solution for A and 5, the determinant must be zero: 



det 



g 4 + ^G lig 2 - u 2 



P ^G 12 q 2 



mi 



V 



P20 , 



o, 



(12) 



/ 



l G 12 q 2 T^q A + m G 22 q 2 - co 2 

which determines the dispersion relation between the frequency and the wave number. One can 
easily check that for the one component condensate case, by putting G\ 2 = 0, Eq.(12) indeed gives 
the usual phonon spectrum. In the present binary BEC mixture case, we have 
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Obviously, both uo\ and uo 2 _ are real as given by Eq.(13). They give the two phonon velocities in 
the binary BEC mixture: 
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However, if 
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the branch uj 2 _ can be negative. This implies an imaginary frequency c<j_ and an imaginary phonon 
velocity, which shows that the initial homogeneously mixed state is unstable. One can verify that 
the sufficient condition for the inequality Eq.(14) to hold for small enough wave number q is the 
validity of the inequality Eq.(4). The modes defined by Eq. (13) with imaginary frequencies will 



then grow exponentially with time. Unlike the usual situation near equilibrium, this growth from 
the present non-equilibrium homogeneously mixed state will be dominated by the fastest growth 
mode. This is precisely the same case as in the initial stage of the classical spinodal decomposition. 

Now we look for the fastest growth mode. It is equivalent to find the most negative value of oo 2 _. 
The condition duj 2 _/d(q 2 ) = leads to 
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Eq.(15) is actually a quartic equation for q^axi whose algebraic solutions can be found by the 
standard way. The corresponding imaginary frequency branch is 
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The detailed general expressions for both q max and cj_, 
ticularly illustrative. We will not explore them here. 

To get a better understanding of the physical implications of Eqs. (15,16), we consider a special 
case relevant to recent experiments where particles of the two condensates have the same mass, 
m i = m 2 = m - In this case we find the wavenumber corresponding to the most negative oj 2 _ max is 
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The physics implied in Eqs. (15,16) or Eqs. (17,18) is as follows. Starting from the initial 
homogeneous mixture of the two condensates, on the time scale given by 

tl 1/1^— max | j 

domain patterns of the phase segregation with the characteristic length 

1 1 1/ Qmax 

will appear. Particularly, for the weakly segregated phase of Gii q 22 — 1 — > 0, we have the length 
scale 



,-1 ^ m 
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and for the strongly segregated phase of Gii q 22 — 1 — ► 00, we have 
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Here Aj = ^/ y Gu/y/GuG^ — 1 and £j = yh 2 /2mj l/pj Gjj are the penetration and healing 
lengths in the binary BEC mixture 0] . Those length and times scales can be measured experimen- 
tally. We will come back to the experimental situation below. 

After the stage I of fast growth into the domain pattern characterized by the length scale Zj, 
the system will gradually approach to the true ground state of the complete phase segregation: 
one condensate in one region and the second condensate in another region. This stage is slow, 
dominated by the slowest mode, and is the subject of next section. 

It is now evident that the stage I of the growth of binary BEC's shares the same phenomenology 
of the initial stage of the classical spinodal decomposition: the domination of the fastest growth 
mode, the appearance of domains of segregated phases, and the conservation of particle numbers. 
There are, however, two important differences. First, the dynamical evolution of the binary BEC's 
is governed by a coupled nonlinear time dependent Schrodinger equations, not by a nonlinear 
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diffusion equation supplemented with the continuity equation, the Cahn-Hilliard equation. There 
is no external relaxation process for the present wave functions. Secondly, the energy of the binary 
BEC's is conserved during the growth process, not as in the case of classical spinodal decomposition 
where the system energy always decreases. 

We wish to point out that the initial stage of spinodal decomposition without diffusion constant 
has been explored in heavy ion collisions for Fermi systems, based on the kinetic equation approach 
@. Qualitatively similar results have been obtained there. The study in alkali BEC mixtures should 
shed light on this important process, because the observation of real time evolution is in principle 
attainable experimentally. 

IV. STAGE II: MERGING AND OSCILLATING BETWEEN DOMAINS 

The BEC binary mixture occurs in a trap. This finite size effect of the droplet leads to the broken 
symmetry of the condensate profiles ||, which tends to separate the condensates in mutually isolated 
regions. This implies that there is no contact between different domains of the same condensates 
formed in the stage I. The classical spinodal process involves diffusion. An estimate of the diffusion 
constant for the BEC system can be made from kinetics theory. The ratio of the time scales for 
quantum and classical particle transport is of the order of the ratio of the BEC cloud size to the de 
Broglie wavelength. This is much smaller than one for the experimental systems of interest. The 
classical diffusion process is thus not important. Because of the domains are not connected and 
because the diffusion for the BEC mixture is extremely low, all the mechanisms for the late stage 
classical spinodal decomposition are not applicable. We propose that it is the Josephson effect that 
is responsible for the approach to equilibrium in stage II. Two models for the Josephson effect, the 
'rigid pendulum' model || and the 'soft pendulum' model |J will be discussed. They both give the 
same time scale when the 'Rabi' frequency is small. Those model studies are performed for a single 
BEC condensate separated by an external potential barrier. We adapt their analysis to the present 
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situation of one condensate sitting in the potential wells formed by another condensate during the 
phase segregation process. 

For a classical spinodal decomposition in a binary fluid mixture, there are several scenarios 
of growth in stage II. (1) For the domains to grow by the diffusion of materials between them, 
a finite diffusion constant is needed. For the binary BEC mixture, the temperature is too low 
for the diffusion to occur. The diffusion constant is practically zero. Furthermore, particles of 
one condensate have to diffuse into the domains of the other condensate because domains of the 
same condensate are not connected. This requires a finite activation energy. Hence this process is 
not possible at low enough temperatures. (2) The same argument also rules out the applicability 
of the diffusion-enhanced collisions by which the diffusion field around the domains leads to an 
attraction between them. (3) Another mechanism is the noise-induced growth. However, there is 
no external noise in the BEC described by the time dependent nonlinear Schrodinger equations. 
This mechanism is thus irrelevant here. (4) A seemingly relevant one is the hydrodynamic growth 
that is driven by the pressure difference between the points of different curvatures. This would be 
a rather fast process. However, this also requires a connected phase of one of the species. Hence, 
at zero temperature and for a finite system, classical scenarios for the later stage will not work 
here. We need to find a scenario working under the present conditions, and at the same time that 
is consistent with the time dependent nonlinear Schrodinger equations. 

The only scenario which can transport particles across a forbidden region at zero temperature is 
the tunneling process. Let us consider the specific case of two domains of condensate 1 separated by 
a domain with width d of condensate 2. The ability of condensate 1 to tunnel through condensate 2 
is described by the penetration depth A, as discussed in Ref. [|J. Hence the probability of condensate 
1 to tunnel through condensate 2 can be estimated as 

p = e~ 2 * , (19) 

when p is much smaller than 1. The finiteness, though small, of the tunneling probability suggests 
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that it is the Josephson effect responsible for the relaxation process in the stage II. The Josephson 
effect leads to the merging of two domains of the same condensate at sufficiently low temperatures. 
The dynamics of such motion may be governed by the 'rigid pendulum' Hamiltonian for a Josephson 
junction ||: 

H{<j>,n) = Ej(l - cos0) + X -E c n 2 , 

where Ej is the Josephson coupling energy determined by the tunneling probability, n = (n x — 
n y )/2 is the particle number difference between the numbers of particles, n x and n y , in the two 
domains, and Ec = dfi/dn is the 'capacitive' energy due to interactions. In the absence of external 
constraints, fi = Ecn. The phase difference between the two domains is conjugated to n, as in 
usual Josephson junctions. Under the appropriate condition, such as low temperature and smallness 
of the capacitive energy, there may be an oscillation between the two domains of condensate 1 
separated by condensate 2. In such a case, we may estimate that the oscillation period tjj = 27i/ujp, 
with the so-called Josephson plasma frequency M 



yEcEj 

uj P = . (20) 

For small tunneling probability, the Josephson junction energy may be estimated as |§ 

Ej = n]^ 3 hu! e~ 2 A , 

and the capacitive energy as 

2 /n T \" - 6 /15a n ^ a4 



e ^-At) Kit) h "<>- 

Here ojq is the harmonic oscillator frequency for condensate 1 in a harmonic trap, ao = Jh/miLUo is 
the corresponding oscillator length, the an is the scattering length of condensate 1, and = n x +n y 
is the total number of particles in domain x and y. Then the oscillatory time scale between the 
domains determines by the Josephson plasma frequency 
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The rigid pendulum Hamiltonian would give a good description when n << t%t- 

Another description of the Josephson effect uses the 'soft pendulum' Hamiltonian proposed in 
Ref. §: 

H((f),n) = -Ej(n) cos (ft + ^E c n 2 , 

with the number dependent Josephson coupling energy Ej{n) = n T Tiuj RS J\ — {n/n T ) 2 . Here lor ~ 
Ej/h is the Rabi frequency determined by the overlap integral between the two wave functions in 
the absence of the mutual repulsive interaction Gu- The appropriate frequency scale in this case is 

u% = u 2 JP + u R , (22) 

which determines the time scale for the two domains to merge. The Josephson plasma frequency 
in Eq.(22) of the soft pendulum model is ojjp = ^tltEcUr/Tl . For the experimental situation of 
interest, ojr « u JP . Thus the two approaches give essentially the same answer. The detailed 
derivation of Ej was given in Ref. ||. 

Given the Josephson effect is the dominant mechanism in the stage II, the time scale to arrive at 
the ground state will be determined by the Josephson effect at the final two domains, in which the 
two domains of condensate 1 is separated by the domain of condensate 2, in the manner of 1 2 1 2 
spatial configuration for the case of equal numbers of the condensates. The width of each domain 
is then D/4, with D the size of the trap. According to the above analysis, the largest time scale 
determined by Eq.(22), the slowest mode, in the stage II is: 

t n = 2n/iu sp . (23) 

The arrangement of 1 2 1 spatial configuration may also occur here, in which it is more likely for 
condensate 2 to tunnel through 1 to the edge of the trap, because of the larger tunneling probability. 
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V. DISCUSSIONS 



The first question is that the quantum spinodal decomposition discussed above can happen 
or not. In terms of the atomic scattering lengths of condensate atoms ajj, the interactions 
are Gjj = Airffajj/rrij. The typical value of ajj for 87 Rb is about 50A. The typical density 
realized for the binary BEC mixture is about p i0 ~ 10 14 /cm 3 . Hence the healing length is 



£ = ^(h 2 /2m)(l/G jj p i0 ) = ^Jl/iSnajjPio) ~ 3000A. For the different hyperfine states of 87 Rb, 
it is known now jl^J^J^I that Gu/ y/GuG^ > 1- Hence the ground state of the phase segregated 
phase can be realizable. Therefore, the quantum spinodal decomposition can happen. 

Experimentally, starting from the initially homogeneous state, after a short period of time a 
domain pattern does appear. Then a damped oscillation between the domain pattern has been 
observed. Eventually the binary BEC mixture sets to the segregated phase flQfl. If we take 



Gu/ VG11G22 = 1.04, the penetration depth is A = £/ 'y Gyij 'VCnG'22 — 1 ~ 1.5/zm. The length 
scale lj determined by Eq.(17) in the stage I is 1.5/zm, which is the same order of magnitude in 
comparison with experimental data |10|. This length is also comparable to the domain wall width 
seen experimentally. The corresponding time scale tj (Eq. (18)) is then 6 ms, again the same order 



of magnitude, though both estimated lj and ti are somewhat smaller. If G12/ \fGuG 2 2 = 1.00 04, 
the corresponding length and time scales in stage I according to the present analysis are lj = 15/mi 
and tj = 600ms, respectively. Time and length scales calculated in this way are larger than the 
numbers extracted from the experiment ||10|| . If we assume the damped oscillation to equilibrium 
observed experimentally is the stage II discussed here, taking the Rabi frequency lor = 0.2nK = 
4 Hz m and the total number of particle Aq = 10 6 , we find the period according to Eq. (23) is 
30 ms, comparable to the experimental value. The estimation from Ref. |J gives a larger value of 
about 200 ms. We think this estimate of the larger time scale is due to the use of a larger estimate 
of d/A. In our view, d and A are comparable in stage I. Thus for stage II, d/A is of the order of 
rid, the number of domains formed in stage I. For the experiment of interest to use, rs 2. Thus 
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for the case at hand, we think the first estimate is more reasonable. 

The present analysis shows that the stage I is insensitive to damping. Similar conclusion has 
been reached in the study of heavy ion collisions for fermions 0. At this moment we do not 
have a reliable estimation of the damping in the Josephson oscillation, as indicated experimentally. 
Nevertheless, given the uncertainty in the value of G\% or a i2 , we conclude that the stage II of the 
quantum spinodal decomposition may have been observed. 



VI. CONCLUSION 

In the present paper we have studied the problem of the dynamical evolution of the binary 
BEC mixtures starting from the homogeneously mixed unstable state. We have found a parallel 
analogy to the usual spinodal decomposition process by calling the present one the quantum spin- 
odal decomposition. The quantum spinodal decomposition consists of two stages: the stage I of 
initial fast growth of domains, characterized by a time and a length scale, and the stage II of a 
slow relaxation towards equilibrium dominated by the Josephson effect. The coupled non-linear 
Schrodinger equations provide a good theoretical description, which enable us to calculate the time 
and length scales. In comparison with recent experiments, we conclude that the quantum spinodal 
decomposition can be, and may have been, realized. 



The length and time scales for Rb mixtures is controlled by the difference of r = G^/v^ii G22 
and 1. Since r is close to 1 experimentally, the data can provide for a very sensitive estimate of r-1. 
Since this parameter determines both the time and length scale of stage I, it is a self-consistent 
check on the physical picture provided here. 

Periodic-like structures have also been observed in the phase segregation of spin 1 Na mixtures 



11]. We think a similar picture of quantum spinodal decomposition applies to that case as well. 
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